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ABSTRACT 

We present the results of an Australia Telescope Compact Array 1.4 GHz spectropolarimetric aper- 
ture synthesis survey of 34 square degrees centred on Centaurus A — NGC 5128. A catalogue of 1005 
extragalactic compact radio sources in the field to a continuum flux density of 3mJy beam" 1 is pro- 
vided along with a table of Faraday rotation measures (RMs) and linear polarised intensities for the 
28% of sources with high signal-to-noise in linear polarisation. We use the ensemble of 281 background 
polarised sources as line-of-sight probes of the structure of the giant radio lobes of Centaurus A. This 
is the first time such a method has been applied to radio galaxy lobes and we explain how it differs 
from the conventional methods that are often complicated by depth and beam depolarisation effects. 
Assuming a magnetic field strength in the lobes of 1.3 B\ fiG, where B\ = 1 is implied by equipartition 
between magnetic fields and relativistic particles, the upper limit we derive on the maximum possible 
difference between the average RM of 121 sources behind Centaurus A and the average RM of the 160 
sources along sightlines outside Centaurus A implies an upper limit on the volume-averaged thermal 
plasma density in the giant radio lobes of (n e ) < 5 x 10~ 5 _Bf 1 cm~ 3 . We use an RM structure func- 
tion analysis and report the detection of a turbulent RM signal, with rms <trm = 17 rad m -2 and scale 
size 0.3°, associated with the southern giant lobe. We cannot verify whether this signal arises from 
turbulent structure throughout the lobe or only in a thin skin (or sheath) around the edge, although 
we favour the latter. The RM signal is modelled as possibly arising from a thin skin with a thermal 
plasma density equivalent to the Centaurus intragroup medium density and a coherent magnetic field 
that reverses its sign on a spatial scale of 20kpc. For a thermal density of n\ 10~ 3 cm~ 3 , the skin 
magnetic field strength is 0.8 n^ 1 /iG. 

Subject headings: galaxies: individual (Centaurus A, NGC 5128) — techniques: interferometric, po- 
larimetric 



1. INTRODUCTION 

The lobes of radio galaxies are magnetised, quasi- 
freely expanding rarified cavities inflated by relativistic 
jets propagating outwards through the intergalactic 
medium from a cent ral supermassive black hole (e.g. 
iBegelman et al.l I1984D . As such, radio lobes could be 
excellent sites for high-energy particle acceleration and 
even the production of ultra-high-energy cosmic rays 
dBenford fc Protherod l2008t iFraschetti fe Melial 120081 : 
Hardcastl e et al.ll2009f) . Knowledge of the physical con- 
ditions in radio lobes, including both the lobe magnetic 
field strength and thermal plasma density, are important 
to e xplore any high-energ y acceleration mechanisms in 
full (|Kronberg et al.ll2004[ ). Understanding the magnetic 
and thermal properties of radio lobes in detail is also 
fundamental to our understanding of galaxy formation 
in terms of feedback processes betwe en the AGN and the 
interstellar / inter galactic med ium (|Croton et alj 120061 : 
ICroft et al.ir2006l : lElbaz et^lWfflh . 

Linearly polarised electromagnetic radiation passing 
through a magnetised thermal plasma causes rotation in 
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the angle of polarisation of the radiation at a rate given 

by, 

69 = 0.81 n e B« SI A 2 

(1) 

= SRM A 2 1 ^ 

where 9 (in radians) is the position angle of the radiation 
at wavelength A (in meters), n e is the thermal electron 
density (in cm -3 ), Bu is the line of sight component of 
the magnetic field (in /iG), I is the path length through 
the rotating material (in pc) and RM is the Faraday 
rotation measure (in units of rad m~ 2 ). 

If radio galaxy lobes contain magnetised, thermal 
plasma they will have an associated intrinsic Faraday 
depth. Observations of RMs in radio sources have shown 
that i nternal Faraday rotat i on in radio lobes is quite 
small dKronberg et al.| | 1986l 120041 : iSchoenmakers et alJ 
1998; Pal ma et al.l 120001 ) with estimates of the thermal 
electron densities of n e <J 10~ 6 cm~ 3 , assuming a 
thermal plasma is distributed uniformly across the 
lobes. Such low inferred thermal matter densities are 
orders of magnitude lower than the upper limits on hot 
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(keV) gas obtained fr om measurements of X-ray cavities 
around radio lobes (IBlanton et all 120011 : IFabian et al.l 
[20001 INulsen et al.ll2005ft ~ 

Whereas only upper limits exist on the uniform 
thermal matter density inside radio lobes, there is 
conflicting evidence regarding the presence of a Faraday 
rotating thin skin (or sheath) around the lobes caused 
by entrainment of interg alactic plasma. For example, 
in the case of Cygnus A, IDreher et all (|1987t ) attribute 
observed RM variations of thousands of rad m~ 2 wholly 
to the foreground i ntracluster med i um th at Cygnus A is 
embedded within. iBicknell et al.l ((1990), however, use 
the same data to show that this RM structure could 
arise from a thin skin around the radio lobes where 
Kevin-Helmholtz instabilities have caused a mixing of 
the lobe plasma with the intergalactic medium. More 
recently, a similar debate has arisen as to the origin 
(intracluster medium versus thin skin) of RMs (and RM 
variations) in excess of ±1000 rad m~ 2 across the radio 
galaxy P KS 1246— 410 in the centre of the Cen t aurus 
cluster 1 (ITavlor et all l2002t iRudnick fc Blundel 1200$ 
lEnsslin et al.ll200a iTavlor et al.ll2007h . 

There are three distinct scenarios we consider (but see 
iBurnl 119661) when using Faraday rotation to probe the 
properties of a magnetised, thermal plasma: 

1. The diffuse polarised synchrotron emission from the 
lobes is mixed with the magnetised, thermal (Faraday 
rotating) plasma. In this case, the emission from 
the back of the source will have been rotated more 
than the emission from the front of the source and 
so 'depth depolarisation' can occur along any line of 
sight. In addition, 'beam depolarisation' can occur 
due to variations in RM on scales smaller than the 
observing beam. In the former situation (i.e. not 
applicable to beam depolarisation) the concept of 
Faraday depth is introduced. Accurate determination 
of Faraday depth is complex, but necessary to extract 
information on the magnetic field and thermal density 
within the source (Ci offi fc Jonesl fl980h . Radio galaxy 
lobes embedded in clusters or grou ps are often used to 
probe the c l uster / g roup medium (IDreher et alj 119871: 
Clarke et all 120011: lEilek fc Owed l2002t ITavlor et al l 
2002HLaing et al.1 12008ft . Here one must first show that 
the rotating plasma arises purely from the foreground 
mediu m itself rather than the lob es or skin of the radio 
source (jRudnick fc Blundellll2003l) . 

2. The polarised emission is diffuse and located 
behind the Faraday rotating plasma. This is similar 
to the above scenario in terms of beam depolarisation, 
however no depth depolarisation occurs because there 
is no mixing of the emitting and rotating regions. For 
example, extended radio galaxies could be used to 
directly probe Galactic magnetic fields. 

3. The polarised emission is unresolved and located 
behind the magnetised, thermal (Faraday rotating) 

1 The Centaurus cluster at z = 0.01 is located behind the Cen- 
taurus group at z = 0.0018 that contains the radio galaxy Centau- 
rus A. 



plasma. In this scenario, all the polarised signal from 
any element of the background source is rotated by the 
entire line of sight through the screen. The screen can 
cause spatial depolarisation due to variations in RM in 
the screen across the angular size of the source. This 
probes scale sizes in the screen on scale sizes smaller than 
the background source size. No depth depolarisation 
occurs. This technique is often used to investigate the 



magnetic structure of the Milky Way dBrown et al 



200S 



20071: iMao et alJl2008ft . nearby galaxies pan et alJH99S 
Gaensler et al.l l2005f) and galaxy clusters (|Kim et al.l 



Note that scenarios 2 and 3 above are not affected 
by back-versus-fr ont differences such as the Lain g- 
Garrington effect (|Lainglll9"8 8: Gar rington et al.l fl988) . 

The typical angular size subtended by radio galaxy 
lobes is too small to include a statistically significant 
number of compact polarised background sources, at 
least for the source densities reached with current 
sensitivity (typically mJy beam -1 ). Hence, up until 
now, all studies of the Faraday rotation in radio galaxy 
lobes have been restricted to using emission from 
the lobes themse lves, as in scenar i o 1 above (recen t 
examples include iKharb et a l. 2009; Lai ng et ail 12008ft . 
The outer lobes of the nearest radio galaxy, Centaurus 
A, subtend a large enough angular size (w 45 deg 2 ) 
that hundreds of polarised sources are detected along 
sightlines behind them. For the first time, we can 
investigate the magnetised plasma in radio lobes - 
using scenario 3 — without the complexities added by 
depth and beam depolarisation effects. This is the basis 
for the analysis presented in this paper. 

We have recently completed a large spectropolarimet- 
ric imaging campaign at 1.4 GHz with the Australia 
Telescope Compact Array (ATCA) and the Parkes 64 m 
radio telescope, to image in full the polarised structure 
of the nearest radio galaxy, Centaurus A. The full spec- 
tropolarimetric images of Centaurus A from ATCA and 
Parkes data combined will be reported in a subsequent 
paper (Feain et al. 2010, in preparation). An additional 
result of the ATCA component of the observations, we 
have also observed in full polarisation 1005 compact 
radio sources in the background of Centaurus A: some 
along lines of sight through the lobes and some along 
lines of sight beyond (outside) the boundaries of the 
radio lobes. In this paper we present a catalogue of 
these 1005 compact radio sources. RMs and polarised 
flux densities are presented for a subset (281) of the 
sources. We investigate the spatial correspondence 
between the radio lobes of Centaurus A and the both 
the distributions of RMs and fractional polarisation of 
the background sources. 

This paper is divided into sections in the following 
way: $2] describes the ATCA radio continuum observa- 
tions and $3] outlines the calibration and data processing 
procedures. In f|4j we define the procedure used for 
source finding, give the format for the source catalogue 
and provide the URL where the entire catalogue can be 
accessed. <J5] describes the Rotation Measure Synthesis 
technique with which we derived reliable RMs and 
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polarised flux densities for 281 out of the 1005 sources 
catalogued. In Sj6] we show the RM distribution and 
briefly compare the total and polarised intensity source 
counts from our data with total and polarised intensity 
source counts from the literature. We also compute 
a lower limit on the very small scale RM fluctuations 
from the lobes of Centaurus A. SJ7] presents a detailed 
investigation into the spatial variations in the ensemble 
of 281 Faraday rotation measures. We fit and subtract 
out the Milky Way foreground RM component leaving a 
residual excess RM dispersion on angular scales 6 ~ 0.3°. 
We model this excess as possibly arising from a thin 
skin around the southern giant radio lobe. Finally, our 
concluding remarks are given <J5] 

We adopt a distance to Centaurus A of 3.8 Mpc from 
iReikubal (|20Q4f l. At 3.8 Mpc, 1° corresponds to «66kpc. 

2. OBSERVATIONS 

The data used in this paper were obtained as part of 
a larger (~45 deg 2 ) radio synthesis imaging survey of 
Centaurus A (NGC 5128), the nearest active galaxy in 
the Universe. 

The ATCA was used in mosaic mode over four 
epochs between 2006 December and 2008 March to 
observe a total field of view covering 5° x 9° (in 406 
pointings) and centred on RA(J2000) 13 h 25 m 27.6 s , 
DEC(J2000) -43°01'09". The standard continuum 
correlator configuration was employed for the obser- 
vations at all epochs. This correlator configuration 
(full_128_2) was chosen because it allows dual- 
frequency observations using 2x128 MHz bandwidth 
observing windows split into 32 channels (the spectral 
resolution is 1.77 channels) per window. We centred the 
two frequency windows on 1344 and 1432 MHz; used 
in this way, the two windows overlap by 40 MHz and 
the total useable bandwith is 192 MHz. Observations 
were carried out with four complementary array con- 
figurations, each with a maximum baseline of 750 m. 
The ATCA's linearly polarised feeds measure cross- 
polarisation data which allowed us to derive the four 
relevant Stokes parameters (I,Q,U,V). The temporal 
variations in the atmospheric phase were tracked with a 
two minute observation of PKS B1316— 46 about once 
an hour. In addition, the parallactic angle coverage of 
PKS B 13 16— 46 was sufficient to enable us to solve for 
the polarisation leakages in each 12 hr observing block. 
PKS B1934— 638 was observed once a day and used 
both to correct for the instrumental bandpass and derive 
the absolute flux density scale. On average, each of the 
406 pointings in the mosaic received 100—120 minutes 
integration. 

3. DATA CALIBRATION AND POST-PROCESSING 

The data were inspected for radio frequency in- 
tereference (RFI) and flagged accordingly using Root 
Mean Square (RMS)-Based flagging in the automate d 
RFI detection algorithm pieflag (|Middelberd 120061) . 
Approximately 5% and 30% of our visibilities at 1384 
and 1432 MHz, respectively, were flagged. Each baseline 
and channel was then visually inspected and, where 



residual RFI was evident, manually flagged, miriad 
was used in the standard way to derive and apply the 
instrumental bandpass, gain, phase and absolute flux 
density calibration. The typical level of polarisation 
leakage in an individual pointing was ~ 0.1% at the 
phase centre, ~ 0.5% at the half-power point and 
~ 1.7% at 0.1-power. The off-axis polarisation leakage 
of the ATCA is reduced substantially by the combi- 
nation of mosaicing and tracking a source over a wide 
range of paralactic angles; we measured the leakage 
across the calibrated final mosaic to be ~ 0.1%. After 
calibration and disregarding the edge channels, there 
were 24x 8 MHz channels between 1296 and 1480MHz. 
Faraday rotation measures were derived across the 
full 24-channcl band. The ionospheric variations in 
RM above the ATCA over the co urse of the observa- 
tions are typically < l rad m ~ 2 IjMacMillan fc Mausl 
120051 : iBilitza fc Reinischl [200l . The variations in the 
ionospheric RM between observations will cause a 
slight amount of depolarisation, leading to a few per- 
cent decrease in the measured polarisation of the sources. 

In this paper, we concentrate on the polarisation 
properties of the compact extragalactic radio sources 
in the background to Centaurus A. For the purposes 
of analysing these background radio sources, the fore- 
ground large scale emission from the giant lobes of 
Centaurus A was considered to be contamination. We 
therefore filtered out this large scale emission from our 
data by disregarding all visibilities from baselines shorter 
than 300 metres (1.4 kA). The final angular resolution 
for this data (after filtering) is 60"x33", 3° (as given 
in Table [J). For this work, each of the 406 ATCA 
pointings were considered, calibrated, deconvolved and 
restored individually (using the same size beam for 
all pointings) and then stitched together with linear 
mosaicing. For each individual pointing, a Stokes /, Q, U 
and V image was produced (the latter to measure the 
noise properties in each channel) for every one of the 
24 independent (8 MHz wide) frequency channels and, 
for Stokes I only, a continuum imag e was also created 
using multi-frequency synthesis (jSault fc Wieringal 
1994). For the latter Stokes I image alone, a single 
iteration of phase-only self-calibration was applied; no 
self-calibration was necessary for the Stokes Q or U 
images. 

A mask was applied to the spatially filtered dataset 
in order to create an image with a roughly uniform 
sensitivity; the point source sensitivities for the final 
continuum image as well as the individual 8 MHz Stokes 
Q and U images are given in Table [TJ Figure [1] shows 
the mask that was applied to the spatially filtered 
image: the large blanked region in the centre of the 
mosaic is due to high residual sidelobes around the 
very bright (S pca k ~ 7.5 Jy beam -1 at 1.4 GHz) central 
region of Centaurus A. The two smaller blanked regions 
south of the Centaurus A core are associated with the 
residual sidelobes of two bright background sources 
(pks 1320-446 and pmn J1318-4620). The edge 
of the mosaic was also masked to remove elevated 
noise levels (;> 0.3 mJy beam -1 ) where primary beam 
correction was significant. Figures [2] and [3] show the 
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TABLE 1 

Observational Parameters 



Parameter 


Continuum 


Stokes (Q, U) 


Bandwidth (MHz) 


192 


8 per channel 


(j (mjy beam -1 ) 


0.15 


0.5 per channel 




63"x33" 


63"x33" 


Position Angle 


3° 


3° 



Position angle is defined with north 0° and east +90°. 



continuum intensity and corresponding linear polarised 
intensity (not corrected for Ricean bias) for a rep- 
resentative portion of the surveyed area. Figure [1] 
has a roughly uniform sensitivity that is given in Ta- 
ble[T] The total survey area, after masking, is 33.93 deg 2 . 

4. SOURCE FINDING 

Source de tection was performe d using the miriad rou- 
tine SFIND (jHopkins et al.1 [20021 on the total intensity 
image in Figure [TJ We ran SFIND in its original mode, 
which uses a 'Search and Destory' algorithm much the 
same as the AIPS task VSAD, to find all sources in the 
image and fit them with a gaussian. This extracted an 
initial list of candidate sources with a peak flux density 
greater than 3mJy beam -1 (20ct). To obtain more accu- 
rate fits than those produced by SFIND we then ran the 
miriad task imfit to do a constrained gaussian fit (re- 
stricting the fit to a small region around the object) for 
each source detected in the initial candidate list. We ran 
this iteratively, subtracting each fitted source from the 
image, to produce a residual map. We identified poor 
fits by comparing the rms noise in a small region around 
the source in the original image and in the residual map. 
Cases in which the rms noise increased after the fitted 
source was subtracted from the image were investigated 
further. Sources with poor fits and sources with multi- 
ple components are identified as such in the catalogue. 
The final catalogue consists of 1005 compact sources to 
a detection threshold of 3mJy beam -1 . The associated 
error in the peak flux density for each source has been 
estimated by the quadrature sum of the rms noise in 
the image and the uncertainty in transfcring the flux 
scale of the ATCA primary calibrator pks B1934-638 
(1— 2% 2 ). The error in the integrated flux density, AJ, 
for each sourc e was estimated, using Equation 16 in 
iCondonl (|1997| ) to be, 

A/ = AAx4, (2) 

where A is the peak flux density. The format of the 
catalogue in Table [2] is: 

Column (1): Source Name. 

Columns (2,3): Right Ascension and Declination in 
J2000 coordinates. 

Columns (4,5): Peak flux density averaged over the 
full bandwidth in units of mjy and its associated error. 
Columns (6,7): Integrated flux density averaged over 
the full bandwidth in units of mjy and its associated 
error. 

2 see ATNF Technical Memo AT/39.3/040. 



Columns (8,9,10): Fitted major and minor axis and 
position angle from imfit 

Columns (11,12,13): Deconvolved major and minor 
axis and position angle from imfit 



5. POLARISATION AND FARADAY ROTATION MEASURE 

The polarised fluxes and Faraday rotation measures 
(RMs) were derived for 281 of the 1005 catalogued 
sources listed in Table [2] as follows. 

At each source position in Table O we extracted the 
Stokes Q and U values and their expected RMS error 
(based on the sensitivity of the map at that position) 
in each of the 24 independent spectral channels in our 
data-set. Converting from frequency to A 2 , where A is 
the observing wavelength of each channel, we then have 
a dependence of the complex Stokes vector (Q, U) as 
a function of A 2 at the peak position of each source. 
For pure Faraday rotation, this vector should have a 
phase that varies linearly with A 2 at a rate equal to the 
source's RM. We extracted this RM from each data-set 
via r o tation measure sy nthesis (|Brentiens fc de Bruvnl 
I2005t iHeald et all 120090 . in which we compute the 
Fourier transform of the complex Stokes vector to yield 
the amplitude of the polarized flux, P, as a function of 
Faraday depth, <j>. For a source with a single valued RM, 
the RM synthesis spectrum will have a single peak whose 
height is equal to the polarized flux, and width (RM 
resolution) determined by the wavelength coverage of 
the observation. The polarisation position angle is given 
by the Q and U values at the peak polarised emission. 
This method achieves signal-to-noise corresponding to 
the full bandwidth of the observation, but with band- 
width smearing effecting only the individual channel 
width. Any deviations from A 2 behaviour resulting from 
complex Faraday rotation structure over the background 
source are seen as structure in the rotation measure 
synthesis spectrum. 

This analysis was applied to all 1005 sources in 
Table [21 with the individual (Q, U) measurements 
weighted by the inverse square of the sensitivity for each 
spectral channel. The total bandwidth and spectral 
resolution of our data mean that the FWHM in Faraday 
depth for a single RM component is 280 rad m~ 2 , 
and that we are sensitive to RMs with magnitudes 
less than «3500 rad m -2 . The resulting Faraday 
depth functions exhibit spectral sidelobes because of 
incomplete wavelength coverage (in the same way that 
an aperture synthesis image shows sidelobes because of 
incomplete u — v coverage). Since we can compute the 
RM transfer function for each source (i.e., the Faraday 
depth spectrum of a source of unit polarized intensity 
and zero RM), we can deconvolve our data-set using 
the same iterative CLEAN app roach routinely applied to 
radio interferometric images (Hogbom 1974). Specifi- 
cally, we have imp lemented RMCLEAN, as described by 
IHeald et all (pool . 

When the above prescription is applied to all sources, 
we now have 1005 deconvolved spectra of P (in units of 
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13M0 m 30 111 20 m 

Right Ascension (JSOOO) 

Fig. 1. — The masked, spatially filtered total intensity radio continuum image of the Centuarus A field used to find and catalogue the 
compact radio sources given in Table [2] The grey shading shows the regions of our field that we used for source finding. In these regions, 
many of the compact radio sources can be seen. Masked regions, shown in white, correspond to areas where the sensitivity of the image 
is poorer, in the vicinity of very bright sources (both the core of Centaurus A as well as two foreground sources with flux densities > 2Jy 
in the southern lobe) or near the edge of the mosaic where primary beam correction was significant. The contours correspond to a Parkes 
1.4 GHz image at 14'resolution (courtesy Mark Calabretta) with levels 1.5, 2, 2.5, 3, 4, 5, 6, 10, 100 Jy beam - 1 . The horizontal-striped region 
beyond the edge of the mask that does not contain point sources is beyond the observed mosaic. The image is shown projected in a SIN 
coordinate system. 
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TABLE 2 

Catalogue of the compact continuum radio sources brighter than 3mJy beam -1 in the field shown in FigureQ] The table 
below shows the format of the catalogue and gives the results for the brightest 20 sources; the full catalogue of 1005 

radio sources is available electronically. 



Source 


RA (J2000) 
hh:mm:ss.ss 


Dec(J2000) 
ddd:mm:ss.s 


Speak 

mjy bm" 1 


ASp Ga k 

mjy bm — 1 


Sint 

mjy 


ASi nt 

mjy 


maj 
II 


mrn 
II 


O 


Tidecon 
maj 
II 


rjdecon 
min 
II 


~D decern 
o 


131452-401056 


13:14:52.67 


-40:10:56.4 


708.7 


7.2 


761.0 


7.7 


63.8 


35.0 


3.1 


11.7 


9.9 


84.8 


132554-464302 


13:25:54.85 


-46:43:02.3 


456.3 


1.6 


545.2 


5.5 


63.4 


39.2 


1.1 


21.3 


6.1 


86.3 


132340-410125 


13:23:40.83 


-41:01:25.9 


441.5 


4.5 


479.1 


4.9 


63.6 


35.5 


3.4 


13.1 


8.5 


82.8 


133130-464508 


13:31:30.39 


-46:45:08.9 


430.0 


1.3 


514.0 


5.2 


63.6 


39.0 


4.0 


21.0 


8.6 


85.8 


131921-443649 


13:19:21.77 


-44:36:49.5 


387.1 


3.9 


453.3 


1.6 


64.1 


38.0 


4.1 


19.2 


11.2 


80.2 


132748-391158 


13:27:48.95 


-39:11:58.1 


355.3 


3.6 


378.1 


3.8 


63.9 


34.6 


2.9 


10.9 


10.3 


-17.1 


132148-383251 


13:21:48.82 


-38:32:51.1 


308.1 


3.1 


328.2 


3.3 


63.8 


34.7 


3.0 


10.7 


10.4 


84.5 


131708-385451 


13:17:08.88 


-38:54:51.4 


302.0 


3.1 


319.0 


3.2 


63.8 


34.4 


3.0 


10.0 


9.8 


-12.0 


132301-384925 


13:23:01.68 


-38:49:25.2 


299.2 


3.6 


321.7 


3.9 


63.8 


35.0 


2.0 


13.1 


8.4 


-52.4 


133118-411956 


13:31:18.68 


-41:19:56.7 


278.6 


1.4 


444.5 


7.1 


65.1 


50.9 


-7.3 


40.2 


12.7 


-75.3 


132104-411451 


13:21:04.51 


-41:14:51.3 


274.4 


2.8 


317.1 


3.3 


66.0 


36.4 


0.1 


22.0 


12.0 


-29.8 


131710-430400 


13:17:10.90 


-43:04:00.5 


267.5 


2.7 


302.5 


3.1 


63.9 


36.8 


3.0 


16.2 


10.9 


-88.0 


132031-410338 


13:20:31.43 


-41:03:38.5 


267.4 


12.9 


326.4 


15.7 


61.8 


41.1 


8.3 








133520-391148 


13:35:20.53 


-39:11:48.4 


253.9 


2.6 


271.7 


2.8 


63.9 


34.8 


3.0 


11.1 


10.7 


-80.1 


132551-400425 


13:25:51.63 


-40:04:25.7 


233.0 


2.4 


250.9 


2.5 


63.7 


35.1 


3.3 


12.2 


9.4 


77.4 


131749-414536 


13:17:49.75 


-41:45:36.1 


226.4 


2.3 


249.9 


2.5 


63.8 


36.0 


3.5 


14.5 


9.8 


80.7 


131739-452221 


13:17:39.58 


-45:22:21.5 


216.0 


2.2 


249.1 


2.6 


63.0 


38.1 


1.1 








132936-442312 


13:29:36.61 


-44:23:12.1 


206.7 


2.1 


237.6 


2.1 


63.6 


37.6 


3.9 


18.2 


8.2 


83.8 


131713-410934 


13:17:13.41 


-41:09:34.5 


197.3 


2.0 


223.6 


2.3 


63.3 


37.2 


3.6 


17.4 


5.6 


86.7 


133357-464138 


13:33:57.33 


-46:41:38.0 


192.9 


2.3 


328.0 


3.9 


88.9 


39.8 


10.8 


63.4 


19.9 


17.1 
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Fig. 2. — A representative portion of the total intensity image 
shown in Figure m The rms level is 0.18 mjy beam . 
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Fig. 3. — As Figure [2] but shown in linear polarised intensity, 
before bias corrrection. The rms level is 0.07 mjy beam -1 

mjy beam ) as a function of <j> (in units of rad m~ 2 ). 
Some examples are shown at varying signal-to-noise in 
Figure LU 

For each source, we identified the peak value of P 
as a function of <fr, and then applied a parabolic fit 
around this peak to yield the best-fit estimate of the 
polarized flux and RM. We then debiased the polarized 
flux by estimating the observed RMS noise, ctqu, in 
the real and imaginary parts of tHe spectrum far from 
the peak, and subtrac ting ctqu in quadrature fr om the 
peak value of P (see iSimmons fc Stewartl I"l985[ ). The 
ratio of the debiased polarized flux to the noise then 
yields the signal-to-noise ratio (SNR) of the detection 
of Faraday rotation. We adopted a threshold SNR > 7 

3 Strictly speaking, P is in units of surface brightness per unit 
</>, i.e., mjy beam -1 rad — 1 m 2 , and the polarized flux is deter- 
mined by integrating P as a function of <j>. However, for sources 
of unresolved Faraday depth the peak surface brightness is equal 
to the total flux, in the same way that the peak surface brightness 
of a point source in an image of intensity is equal to its integrated 
flux. 



as a minimum criteri on for a reliable RM determin ation 
(see, e.g., Figure 9 of lBrentiens fc de Bruvnl feOOS). For 
sources with polarized fluxes above this threshold, we 
computed the uncertainty in RM as the FWHM of the 
RM transfer function divided by twice the SNR. 

Of the 1005 catalogued continuum radio sources, 
281 (28%) polarised sources were robustly detected 
according to the detection criteria described above and 
are listed in Table [3l The format of Table [3] is: 

Column (1): Source Name. 

Columns (2,3): Right Ascension and Declination in 
J2000 coordinates. 

Columns (4,5): Peak polarised flux density (after 
correction for the Ricean bias) in units of mJybm -1 
and its associated signal to noise ratio (SNR). 
Columns (6,7): Measured Faraday rotation measure 
and its associated uncertainty derived using RM synthe- 
sis and RM clean in units of rad m -2 . 
Columns (8): Column 6 after correcting for the 
Galactic contribution (see [JFJ in units of rad m -2 . 

Figured] shows the distribution in the total intensity to 
polarised intensity plane of the 281 sources in our sam- 
ple whose debiased polarised intensities have a SNR>7. 
The dashed, diagonal lines represent lines of constant 
fractional polarisation. 

5.1. Which sources are behind the lobes? 

For much of the rest of this paper, we wish to anal- 
yse separately the polarised sources located behind the 
lobes of Centaurus A from polarised sources whose sight- 
lines pass outside of the lobes. The latter are used as a 
control sample and to model and subtract the RM com- 
ponent from the Milky Way. We define what consti- 
tutes behind the radio lobes to correspond to an 'edge' 
where, at 1.4 GHz and 14' resolution, the surface bright- 
ness of the lobes drops below 1.5 Jy beam -1 ; this is 
the last contour in Figures [T] and [7l At approximately 
this value, the surface brightness of the lobes decreases 
sharply (but there is clearly still diffuse emission from the 
lobes present, probably out to the edges of our mosaic). 
Using a 1.5 Jy beam -1 threshold, 121 out of the total 281 
RMs are behind the radio lobes with the remaining 160 
RMs along sightlines outside the lobes. We have tested 
the sensitivity of our results to the exact value of the 
boundary chosen for inside versus outside the lobes. We 
are confident that our result is robust for all contour val- 
ues that we tested as boundaries (which were between 1 
and 2 Jybeam -1 ). 

6. RM DISTRIBUTION & SOURCE COUNTS 

6.1. RM Distribution 

A histogram of the distribution of the full RM sample 
is shown in the top panel of Figure [6] The mean RM of 
the total distribution is —57 rad m~ 2 with a standard 
deviation of 30 rad m~ 2 . These values are in good 
agreement with earlier observations of the Faraday 
rotation measure of the lobe emission of Centaurus A 
reporte d to vary from —40 to —70 rad m~ 2 a cross the 
source (jCooper et aIJll965t iClarke et al.lll992t ). Other 
studies have also reported an average RM for Centau- 
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Fig. 4. — Results of rotation measure synthesis for eight of the 1005 radio sources in Table [2] shown in increasing order of signal-to-noise 
ratio (SNR) in linear polarised flux density. Each panel plots polarized flux, P, vs. Faraday depth, <j>, at the peak pixel of a source. The 
solid curve shows the amplitude of the complex polarization vector, while the dotted line shows the 1-cr RMS noise level at that position. 
The legend to each panel provides the source name and the SNR, along with the debiased polarized flux and rotation measure in cases 
where SNR> 7. 



TABLE 3 

Catalogue of the 281 polarised radio sources detected. The columns in this table are described in sJU The table below 
gives the brightest 20 sources in linear polarised intensity; the full list of 281 radio sources is available electronically. 



Source RA (J2000) DEC (J2000) P SNR RM ARM* RM** r 







hh:mm:ss.ss 


ddd:mm:ss.s 


(mjybm x ) 




rad m 2 


radm 2 


radm 2 


133357- 


-464138 


13:33:57.34 


-46:41:38.04 


13.95 


154.2 


-59.5 


0.9 


1.6 


131713- 


-410934 


13:17:13.42 


-41:09:34.56 


12.98 


123.4 


-74.6 


1.1 


-19.6 


132937- 


-444416 


13:29:37.75 


-44:44:16.44 


12.47 


129.9 


-77.7 


1.1 


-20.8 


133628- 


-403529 


13:36:28.46 


-40:35:29.40 


10.06 


95.9 


-15.2 


1.4 


20.6 


131749- 


-414536 


13:17:49.75 


-41:45:36.00 


9.43 


76.0 


-30.5 


1.8 


25.9 


132543- 


-461919 


13:25:43.01 


-46:19:19.56 


9.14 


122.9 


-47.2 


1.1 


19.0 


133158- 


-410203 


13:31:58.13 


-41:02:03.12 


8.77 


108.0 


-66.6 


1.3 


-25.3 


133027- 


-452406 


13:30:27.55 


-45:24:06.12 


8.06 


84.7 


-62.3 


1.6 


-3.6 


132553- 


-462121 


13:25:53.71 


-46:21:21.96 


7.79 


94.7 


-42.2 


1.5 


24.0 


133800- 


-452717 


13:38:00.26 


-45:27:17.64 


7.66 


116.6 


-40.7 


1.2 


12.3 


131921- 


-443649 


13:19:21.77 


-44:36:49.68 


7.52 


67.4 


-73.0 


2.1 


-24.7 


133241- 


-450802 


13:32:41.71 


-45:08:02.04 


7.52 


101.9 


-80.6 


1.4 


-7.9 


132145- 


-403618 


13:21:45.46 


-40:36:18.36 


6.53 


90.1 


-71.2 


1.5 


-22.4 


133001- 


-413824 


13:30:01.15 


-41:38:24.00 


6.42 


51.5 


-53.2 


2.7 


-8.0 


132827- 


-464808 


13:28:27.74 


-46:48:08.28 


6.20 


89.2 


-47.7 


1.6 


18.1 


132253- 


-461330 


13:22:53.52 


-46:13:30.72 


6.12 


71.9 


-57.7 


1.9 


10.4 


133746- 


-451641 


13:37:46.54 


-45:16:41.52 


6.05 


69.6 


-56.8 


2.0 


-4.3 


131350- 


-410134 


13:13:50.95 


-41:01:34.32 


6.02 


29.6 


-67.0 


4.7 


-9.4 


131630- 


-444944 


13:16:30.31 


-44:49:44.76 


5.79 


72.8 


-97.9 


1.9 


-29.5 


133107- 


-465745 


13:31:07.06 


-46:57:45.36 


5.57 


69.8 


-18.9 


2.0 


45.4 



* The associated uncertainty in RM does not include the uncertainty due to ionospheric variations which are < 1 rad m 1 . 

** RMs after correcting for the Galactic contribution. 



rus A of -59 ± 3rad m" 2 dGardner fc Daviesl 119661 : 
iSimard-Normandin et~aH 1 1 98 lf > . The RM of Centaurus 
A is typical of that measured for other nearby sources 
which implies that the RM is dominated by a foreground 
comp onent (|Johnston-Hollitt et alj l2004t iShort et al.1 
l2007h . The middle and bottom panels of Figure |6] show 
histograms of the RM distributions of sources behind 
and outside the lobes of Centaurus A, as defined in i l5.ll 



Examination of the available RM sky data from inter- 
polat ed maps (| Johnston-Hollitt et al.1 120041 iFrick et ail 
2001) shows the that Centaurus A is embedded in a large 
region of negative rotation measure of average value ~ 
-27 rad m~ 2 with a slope of roughly -0.9 rad m~ 2 per 
degree in longitude and -0.6 rad m~ 2 per degree in lati- 
tude over the region corresponding to these observations. 
Unfortunately, the paucity of rotation measure data in 
the southern hemisphere means interpolated values in 
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Fig. 5. — Total versus linearly polarised flux density plane for 
the 281 sources in our sample whose debiased polarised intensities 
correspond to SNR>7; see j5]for a justification of this threshold. 
The dashed, diagonal lines represent lines of constant fractional 
polarisation (from top left to bottom right: 100%, 10%, 1%, 0.1%). 

this region have been inferred from some of the the most 
sparse data distributions in the entire sky. 

6.2. Source counts 

We have compared our total and linearly polarised in- 
tensity source counts as a function of flux to published 
source counts in order to confirm that our sample is sta- 
tistically comparable to other samples in the literature 
and that Centaurus A in the foreground is not affecting 
the global statistics of the sample. We find: 

1. Our derived total intensity source counts are con- 
sistent, within the errors, with the continuum 
source counts from the NRAO VLA Sky Survey 
(|Condon et al.lll99l . 

2. Our polarised intensity source counts (in the flux 
range 0.5 < P < 14mJy) are consistent within 
the errors with the polarised source counts from 
the ELAIS Nl field (in the flux range .5 < P < 
20mJy) reported bv lTavlor et"aH (|2007f ). 

6.3. A limit on small scale turbulence in the lobes 

Random fluctuations in either B\\ or in n e on some 
angular scale R will increase the dispersion in RM 
determined on this same scale. If R is smaller than 
the resolution element used to probe the fluctuations, 
then this effect will lead to 'beam' depolarisation where 
the 'beam' size is the angular scale of the background 
polarised sources. Such small scale RM fl uctuations have 
been inferred in the Galactic plane (lHaverkorn et alJ 
2008) and the Large Magellanic Cloud (jGaensler et alJ 
20051 ) using an ensemble of extragalactic polarised 
sources. 

We find no significant difference in the polarised source 
counts for the 121 sources behind the radio lobes of Cen- 
taurus A compared with the 160 sources outside the 
lobes. The lack of detection corresponds to a limit on 
magnetised thermal structure in the lobes of Centaurus 
A on scales much smaller than the median angular scale 
size of the background sources t hemselves; this is i ~ 10 " 
for a lOmJy source at 1.4 GHz (| Windhorst et al.lll984T ). 




-100 -50 
Ratal, on Measure (rad/m"2) 

Fig. 6. — Top panel: RM distribution of the 281 sources in Ta- 
ble [3] The mean RM is —57 radm -2 and the standard devia- 
tion is 30 radm" 2 . Middle panel: RM distribution of the 121 
sources behind the lobes (mean —61.3 radm -2 and standard de- 
viation 29.8 radm -2 ). Bottom panel: RM distribution of the 160 
sources outside the lobe (mean —52.9 radm -2 and standard devi- 
ation 29.2 radm -2 ). The bin size is 20 radm -2 in all panels. 

The mean fractional polarisation of the 121 polarised 4 
sources behind the lobes (quoted with the standard er- 
ror in the mean) is 7.6 ± 0.5% compared to the mean 
fractional polarisation f or the source s outside the lobes 
of 7.3 ± 0.4%. Using the|Rurn| (jl966f ) law for depolarisa- 
tion, 



fd 



(3) 



and a conservative limit of fd > 0.7 on the ratio of the 
observed to intrinsic polarisation of sources behind the 
lobes, we derive an upper limit on the very small scale 
RM turbulence in the giant radio lobes of Centaurus A 
of < 10radm~ 2 on scales <Cl80pc. A similar limit holds 
for the volume-averaged RM amplitude in the radio lobes 
f ^7.ip . implying that within the limits of our data the 
ordered component of the lobes has no prefered scale 



7. SPATIAL VARIATIONS IN THE ROTATION MEASURE 

The spatial distribution of background RMs is repre- 
sented in Figure [71 The size of the disks correspond to 
the value of the RM relative to the mean of the sample 
(— 57rad m~ 2 ). Black and white sources are those 
with positive and negative residuals from the mean, 
respectively. 

A non-zero RM component associated with the lobes 
of Centaurus A could result from thermal material, 
mixed with non-thermal material, in the lobes or, 
perhaps, from the thin skins of ionised plasma around 
the edges of the radio lobes that are created when the 
radio lobes inflate cocoons, sweeping up the thermal 
plasma as they expand into the intergalactic medium. 
There will be an excess RM if the plasma retains an 
ordered component and an excess RM dispersion in the 

4 Note that we analyse those polarised sources with a SNR>7 
only; there are polarised sources below this SNR cutoff but we do 
not consider them in this study. 
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13 h 40 m 30 m 20 m 10 m 

Right Ascension (J2000) 

Fig. 7. — Locations and RMs of the 281 sources in Table [3] To better highlight the variations, the diameter of the sources represent 
the amplitude of their residual RM after the mean RM of the whole distribution (— 57rad m -2 ) has been subtracted. Black and white 
sources are those with positive and negative residuals from the mean, respectively. Overlaid are Parkes 1.4 GHz radio continuum contours 
of Centaurus A. Contour levels are 1.5, 2, 3, 4, 5, 6, 10, 100 Jy beam -1 . The legend on the right hand side of the figure shows the relation 
between the source diameter and the absolute value of the mean-subtracted RM in units of rad m — 2 . 



Faraday rotation structure in the lobes of Centaurus A 



11 



presence of a turbulent component. 

In this section, we are interested in separating the 
Milky Way RM component from any RM component 
that is intrinsic to the radio lobes of Centaurus A. Al- 
though it is not possible to do this for the RM of an 
individual source, it is possible to isolate the Milky Way 
and Centaurus A RM contributions in a statistical sense, 
subject to the assumption that the statistics of the Milky 
Way spatial RM fluctuations on small scales (less than a 
few degrees) do not change appreciably across the field. 
To do this we use the 160 RMs that are not located along 
sightlines through the radio lobes to fit for the fore- 

ground Milky Way component. This was achieved using 
the miriad task impoly to fit a first order polynomial to 
the RM surface in two dimensions. This surface was then 
subtracted from all 281 RMs. In this way, we have essen- 
tially also removed the Milky Way RM component from 
the RMs located behind the radio lobes of Centaurus A. 

7.1. The ordered component 

The measured RMs of the 121 sources behind the 
lobes of Centaurus A have a mean of —61.3 radm -2 
and a standard deviation of 29.8 radm -2 . After re- 
moving the Milky Way RM component using the pro- 
cedure described above, the mean (quoted with the 
standard error in the mean) of the sources reduces to 
— 4.5 ± 2.8 radm" 2 and the standard deviation increases 
slightly to 30.4 radm -2 . The measured RMs of the 160 
sources located along sightlines outside the lobes have 
a mean of —52.9 radm -2 and a standard deviation of 
29.2 radm -2 . Following subtraction of the Milky Way 
RM component, the mean reduces to — 1.6±2.1 rad m -2 
and the standard deviation decreases to 26.6 radm -2 . 
The 3<T upper limit on the difference between the mean 
RM inside and outside the lobes (~ 10 radm -2 ) is con- 
sistent with the limit infered for very small scale RM 
structure in the lobes f "i6.3p . implying that within the 
limits of our data the ordered component of the lobes 
has no prefered scale size. To estimate the uniform ther- 
mal plasma density in the lobes, we assume that the 
lobe magnetic field is dominated by an uniform com- 
ponent with strength B» = 1.3 B\ /xG (setting B\ = 1 
gives the equipartition field strength; lHardcastle et al.1 
2009) and N possible reversals of the field, the path 
length through the lobes to be D = 200 kpc, and the 
upper limit on the volume-averaged RM amplitude as- 
sociated with the lobes to be 10 radm -2 . Substituting 
these values into Equation [TJ we derive an upper limit 
on the uniform thermal plasma density in the lobes of 
(n e ) < 5 x lO -5 .^ -1 \/N cm -3 . For N = 1, this limit 
is 2 — 3 times smaller than previously publ ished lim- 
its (jMarshall k. Clarld Il98ll: llsobe et all 120011) which in 
turn further constrains the expected 7-ray emission in 
the lobes and the prospe cts for the production of ultra- 
high energy cosmic rays (jHardcastle et al.| [2009) . 

7.2. The turbulent component 

The difference between the standard deviation of the 
foreground-subtracted RMs behind the lobes compared 
to outside the lobes (3.8 rad m -2 ), if real, implies struc- 
ture/turbulent magnetised plasma in the lobes. To verify 
and quantify this difference, we have compared the am- 
plitude of the RM fluctuations as a function of angular 



separation in the lobes (or around the edges of the lobes) 
to a 'control group' of RMs located on sightlines outside 
the lobes using an approach based on the structure func- 
tion, 

D R u(e) = ([RM(0 + 0')-RM(0')] 2 ). (4) 

where is the angular separation between sources. The 
structure function is a robust and reliabl e means of mea- 
surin g rotation measure fluctuations iSimonetti et al.l 
SHI . In comparison with other statistical techniques 
like an autocorrelation function analysis for example, the 
structure functions is far less susceptible to uncertainties 
in the mean RM level and to large-scale RM gradients, 
on scales below those comparable to the size of the 
gradient. A structure function-based analysis is also 
immune to the irregular spatial sampling of our RMs. 



7.2.1. Rotation measure structure functions 

Prior to forming structure functions for sources behind 
and outside the lobes, the large-scale RM component 
from the foreground Milky Way was fit and subtracted, 
as described above. When we compared the structure 
function formed from the foreground-corrected RMs 
to the structure function formed from the observed 
RMs (i.e. with no foreground subtraction applied), we 
found no significant difference between the two results 
for scales below 4°. This is unsurprising because 
the large scale gradients are expected to affect the 
structure function on scales comparable to the scale of 
the gradient only. 

Rotation measure differences, as probed by a structure 
function, can arise from intrinsic differences in the 
radio sources themselves, or due to spatial variations 
in the foreground (i.e. in the Galaxy or in the lobes 
of Centaurus A). Figures and [8b are the rotation 
measure structure functions for sources located behind 
and outside the radio lobes of Centaurus A, respectively. 
It is clear from a comparison of Figures [8^, and [8}d that 
the RMs behind the lobes have excess structure on scales 
0.2° <J 9 <J 4.0°. The point at the smallest angular scale, 
9 = 0.025°, is determined solely on the basis of rotation 
measure differences between each of the double sources 
in our sample. Since many of these double sources are 
likely to be the lobes of an individual galaxy, variatons 
in RMs between the double sources is likely to be much 
lower than intrinsic variations between independent 
sources. Thus the jump in amplitude between angular 
separations 9 = 0.025 and > 0.075° is partially due to 
the fact that the latter set is necessarily determined 
using independent sources, with independent intrinsic 
RMs. The variation in the intrinsic (internal) RMs 
of independent sources contributes a "white-noise" 
component to the RM variability, which manifests it- 
self as an offset in the amplitude at all points 9 > 0.075°. 

From the definition of the structure function in 
Equation 01 it can be easily shown that the difference of 
two structure functions is equal to the structure function 
corresponding to the difference of two equivalent power 
spectra. Furthermore, since the intrinsic RM variations 
of the sources are independent of any foreground struc- 
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ture, taking the difference between the two structure 
functions removes the offset. As such, subtracting the 
RM structure function outside the lobes from the RM 
structure function behind the lobes (i.e. Figure — 
Figure [5Jd) gives us a direct measure of the structure 
function of RM fluctuations due to Centaurus A alone. 
This is shown in Figure |Sfc. 

Figure[8]; shows a significant signal at 9 > 0.2°, which 
implies that the lobes of Centaurus A contribute a mea- 
surable RM structure only on scales greater than this. 
The signal saturates at 9 m 0.3° and, within the present 
errors, is consistent with a constant amplitude thereafter. 

We performed a number of additional tests on the data 
as a means of establishing the reality of this signal and 
isolating its origin. The most obvious cause of a spurious 
signal is that there is some large scale asymmetry in 
the spatial rotation measure fluctuations between the 
northern and southern lobes due to a gradient in the 
properties of the Milky Way. Figures [9^-c explore this 
possibility and confirms that the signal in Figure [Sfc 
arises due to RM fluctuations associated with the 
southern lobe of Centaurus A. Figure [SJa compares the 
RM structure function of all the polarised sources that 
are located south of declination —43° to all the polarised 
sources located north of declination —43°. There is 
a clear strong north-south asymmetry present in the 
data, indicating an excess of RM structure in sources 
in the south relative to those in the north. Figure^ 
compares all sources that are outside the lobes and 
located south of 8 = —43° to all sources located outside 
the lobes and north of 8 = —43°. Figure[9j3 shows no 
evidence for a north-south asymmetry outside the lobes. 
FigurelHt compares sources located inside the southern 
lobe to those inside the northern lobe. Here again 
the asymmetry between the lobes is apparent but the 
significance is degraded by the smaller number of sources 
within each lobe. Taken together, Figures [9^-c suggest 
that the north-south asymmetry is due to variations 
between the northern and southern lobes of Centaurus 
A, and not due to the Milky Way. 

Note that the source density in each region is accounted 
for in determining the error bars associated with that 
structure function and does not otherwise systematically 
bias the structure function. 

7.3. Physical interpretation of the lobe RM 

The total and polarised intensity structure of Centau- 
rus A correlates very well with the RM variations of the 
background source s (Feain et al. 2010, in preparation; 
Uunkes et alll993D . Such a correlation could not occur if 
the RM s tructure is un r elated to Centaurus A itself. For 
example, Uunkes et al.l (|1993l ) investigated the polarised 
structure of Centaurus A at 6.3 cm and showed there is 
a striking difference in the polarised intensity structure 
between the northern and southern lobes. The polarised 
emission in northern lobe largely follows the continuum 
emission uniformly down to the sensitivity limits of the 
survey. There is little evidence for depolarisation in 
the northern lobe. The e mission in the southern lobe , 
however, is depolarised (|Gardner fe Whiteoakl 119711) 
and chaotic with position angle jumps of up to 90° in 



places. Between 5 = -44°45' and 8 = -45°45' (J2000 
coordinates) the southern lobe becomes highly c haotic 
and turbulent (see Figure 3b in Uunkes et al.| [1993). The 
good spatial coincidence between the depolarization 
of the southern lobe and that of the RM signal in 
the background sources on angular scales between 
0.2° 9 <^ 4.0° inside the lobes gives us confidence that 
the depolarising (rotating) medium is intrinsic to the 
southern radio lobe of Centaurus A. 

The amplitude of the RM signal in Figure [Sfc is 
approximately 620 rad 2 m~ 4 , implying that the lobes 
of Centaurus A contribute an RM signal with an rms 
gjim ~ 17 radm~ 2 on an angular scale ^0.3° (20kpc 
at the distance to Centaurus A). This RM signal could 
arise either from a thin skin around the lobe or a 
turbulent medium throughout the lobe. Our data are 
not sufficient to confidently distinguish between these 
two possibilities, but we tend to favour the former (thin 
skin) based on the location of the maximum dispersion 
of RMs seen in Figure [7| being well aligned with total 
intensity features along the western boundary of the 
southern lobe (Feain et al. 2010, in preparation) which 
look like surface-wave instabilities. We model this signal 
as follows: 

Suppose that the edge of the southern lobe has uniform 
thermal density, n^ m , and a coherent magnetic field, 
B tl with angle of orientation, 9 ri which changes direction 
randomly on a s patial scale of / ~ 20 kpc. Then, following 
IGaensler et al.l (|2001f ). the dispersion of the RMs pro- 
duced by unformly polarised emission passing through 
the edge of the southern lobe is given by: 



°~ RM 
rad m _ 




where D = 180 kpc is the estimated path length through 
the edge of the lobe. To estimate n^ m we assume 
mixing occurs between the lobe magnetic field and 
the Centaurus intra-gr oup medium due t o Kevin- 
Helmholtz instabilities (jBicknell et alj [1990). Based 
on the ram pressure stripping arguments to explain HI 
gas depletion in dwarf galaxies in the Centaurus group, 
IBouchard et al.l (|2007f ) infer an intragroup medium den- 
sity near Centaurus A of n e ss 10 ~ 3 cm~ 3 ; note however 
that based on X-ray observations. [Feigelson et alj (|198lD 
estimate a similar density for the interstellar medium 
of NGC5128 itself out to 9 kpc. Substituting a RM ,D 
and I into Equation [5] results in Bt ~ 0.8 n^ 1 fiG, where 



-. skin 
l th 



m io- 



The possible existence of a magneto-ionic skin around 
the edges of Centau rus A — a low power F anarof-Rilcy 
class I (FRI) source (|Fanaroff fc Rilevl[l97l — with RM 
fluctuations on scales of 20 kpc, is not dissimilar to that 
inferred for skins around the lobes of the powerful FRII 
sources PKS 2104-2 5N and Cygnus A (jBicknell et all 
119901: ICameronlH988l ). In the latter two cases, the RM 
signals are between 10 and 100 times larger than that 
seen in Centaurus A, which is not unexpected given both 
sources are located at or near the centres of clusters 
where the intracluster density will be higher, possibly 
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Fig. 8. — The RM structure functions of (a) the 121 polarised sources located behind the radio lobes of Centaurus A and (b) the 160 
polarised sources located along sightlines outside the radio lobes of Centaurus A. (c) The difference between (a) and (b). 
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Fig. 9. — (a) The RM structure function of the sources that are located north of S = —43° subtracted from the RM structure function 
of all sources that are located south of <5 = —43°. (b) The structure function of all sources outside the lobes that are located north of 
5 = —43° subtracted from the RM structure function of all sources outside the lobes that are located south of S = —43°. (c) The RM 
structure function of all sources behind the northern lobe subtracted from the RM structure function of all sources behind the southern 
lobe. 



by several orders of magnitude. 

8. CONCLUSION 

In this paper, we have presented the results of a 
1.4 GHz spectropolarimetric imaging survey of 34 deg 2 
centred on the nearest radio galaxy Centaurus A. A 
catalogue of the 1005 background, compact radio sources 
brighter than 3mJy bm has been made available in 
electronic format. 

We used Faraday rotation measure (RM) synthesis to 
measure linear polarised intensities and RMs for each 
source in our sample. There were 281 sources, out of 
1005 in total, with a SNR>7 in linear polarised intensity 
and for these we also publish a table of polarised inten- 
sities and RMs. The density of RMs published here, 
8.3 per deg 2 , represents the densest RM grid known to 
date. Of the 281 RMs, 121 come from sources located 
behind the radio lobes of Centaurus A; the remaining 
160 being from sources located along sightlines outside 
(in projection) the radio lobes. 

The continuum source counts and the linear polarised 
intensity (P) source counts (for sources with SNR>7 in 
P) are consistent with published source counts giving 
us confidence that our sample is statistically robust and 
the emission from the lobes of Centaurus A has not 
contaminated the global statistics of our sample. There 
is no measurable difference between the polarised source 
counts behind the lobes and those outside the lobes. We 
derived an upper limit on the very small scale turbulence 
in the lobes of < 10radm~ 2 on scales <;180pc. This is 



similar to the limit derived on the volume-averaged RM 
signal from the lobes. 

We modelled the Milky Way contribution to the RMs 
by fitting a first order polynomial surface to the 160 
RMs along sightlines outside the lobes of Centaurus A 
and subtracting this surface from all 281 RMs in our 
sample. The residual RMs behind the lobes have a 
mean of —4.5 ± 2.8 radm -2 and a standard deviation of 
30.4 radm -2 , compared to a mean of —1.6 ±2.1 rad m -2 
and standard deviation of 26.6 radm -2 for RMs outside 
the lobes. Assuming a magnetic field strength in 
the lobes of 1.3-Bi/zG, we derived an upper limit to 
the uniform thermal plasma density in the lobes of 
(n e ) < 5 x 10~ 5 B\ cm~ 3 . This is a factor of 2 — 3 times 
smaller than previously published limits and further 
constrains the prospects for high-energy acceleration 
processes in the lobes of Centaurus A. 

We used a structure function analysis to measure an 
excess RM dispersion in the southern lobe of Centaurus 
A of <trm ps 17 radm -2 on angular scales of ~ 0.3° 
associated with the southern lobe of Centuarus A. The 
RM signal is modelled as possibly arising from a thin 
skin with a thermal plasma density equivalent to the 
Centaurus intragroup medium density and a coherent 
magnetic field that reverses its sign on a spatial scale of 
20kpc. If the Centuarus intragroup density is of order 
10~ 3 cm~ 3 , the skin magnetic field strength is of order 
1 /xG. If n e is an order of magnitude smaller, B t will be 
10 x larger. 
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A more sensitive survey of the southern lobe is now 
required to verify if the RM signal is associated with the 
entire southern lobe, or a thin skin (sheath) around the 
edges, although the latter seems more likely. In addition, 
a deep X-ray observation of the edge of the southern lobe 
would be incredibly revealing in terms of the detailed 
physics it could probe. The southern lobe of Centaurus 
A has typically been neglected in favour of investigations 
of the northern lobe; this is almost certainly because of 
the very interesting northern middle lobe (NML) region 
(there has been no southern counterpart to the NML 
detected). In this paper, we have shown that the edge 
of the giant southern lobe, rather than any part of the 
giant northern lobe, is probably associated with fluctu- 
ating RMs on scales of 20kpc. It is extremely timely 
now to follow-up this region with a multiwavelcngth 
campaign at, in particular, X-ray and deep narrow- 
band optical wavelength regimes. We are currently 
completing a detailed analysis of the full polarisation 
radio continuum images of Centaurus at a resolution 
of « 40". In this analysis, we will be able to compare 
the fluctuating RM region with the polarised and total 
intensity structure of the lobes at a factor of about 5 
times the resolution that has been previously possible. 
If our model of a thin skin is correct, these new images 
should reveal the signatures of surface waves around 



the lobes, which would be evidence of mixing between 
the lobe emission and the Centaurus intragroup medium. 
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